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Background: The skeletal muscle echo intensity (EI) during ultrasound imaging has been 
investigated to evaluate the muscle quality. However, EI fluctuates according to the scanning 
conditions. 

Methods: The motor functions and ultrasound images of 19 elderly (73±3.2 years) and 
19 young (22+1.5 years) individuals were investigated and an EI frequency component was 
assessed for more reliable evaluations. Healthy elderly and young subjects participated in this 
study. The motor functions were assessed during walking and according to the knee extension 
muscle strength. The muscle thicknesses of rectus femoris (RF), vastus intermedius (VI), and 
quadriceps femoris (QF) were investigated. Els were calculated and the mean frequencies of 
the regions of interest (MFROIs) for RF and VI were analyzed. 

Results: Els and MFROIs were greater in elderly subjects than in young subjects (P<0.01 
for RF, and P<0.001 for VI, in Els; and P<0.01 for RF, and P<0.05 for VI, in MFROIs). In 
young subjects, EI of RF was greater than that of VI; however, there was no difference between 
the RF and VI MFROIs in both elderly and young subjects. Els of VI exhibited a significantly 
negative correlation with the QF thickness in both elderly and young subjects. RF MFROIs 
negatively correlated with the QF thickness and positively correlated with EI of VI in elderly 
subjects alone. 

Conclusion: These findings suggest that MFROIs of elderly individuals would have a larger 
value than those of young individuals; moreover, MFROIs did not fluctuate greatly with the tissue 
depth and scanning conditions. MFROIs might be thus useful for further investigations of muscle 
quality and applications for the early prevention of age-related motor functional decline. 
Keywords: elderly individual, motor function, muscle strength, ultrasound, echo intensity, 
frequency analysis 

Introduction 

Motor function is required for most activities of daily living, including locomotion. 
The initial stages of motor function decline are exacerbated in elderly individuals in 
whom these early declines lead to further muscular disuse. This age-related decline 
includes decreased muscle mass (eg, sarcopenia) and can accompany the accumulation 
of connective tissues (eg, adipose) in the affected muscles. 12 The muscle mass can be 
measured using imaging methods such as magnetic resonance imaging or computed 
tomography. 

Given its safety and accessibility, ultrasound has become a widely used imaging 
modality in many medical applications. Decreased muscle thickness indicates decreased 
muscle mass, and increased echo intensity (EI) of the muscle indicates increased 
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adipose level. 3-5 Cadore et al reported a relationship between 
EI and cardiovascular performance. 6 Computer-aided gray- 
scale ultrasound images have been analyzed to evaluate 
the skeletal muscle quality. However, EI fluctuates with 
changes in the applied location or angle of the ultrasonic 
probe relative to the observed organ. Moreover, ultrasound 
devices can modify Els at different values to improve image 
quality or resolution. In ultrasound images, deeper tissues are 
usually modified using a larger EI value in order to adjust 
the echo attenuation. For these reasons, unstable EI values 
for observed tissue quality estimations are very difficult to 
standardize. 

The frequency analysis method has been used to 
evaluate the value qualities for different fields of study. 
For example, in electromyography, frequency analysis 
has been used to evaluate the muscle activity quality. 7 - 8 
Frequency analysis of Els of some muscle tissues might 
represent a more reliable method to evaluate the quali- 
ties of such tissues. Fast Fourier Transform (FFT) is the 
most general frequency analysis method; however, it was 
designed on the assumption that the calculated data range 
from "infinity to infinity. Before analysis, the data are 
usually modified using some window function to correct 
the gap between the imaginary infinite data and real finite 
data if the amount of data is not sufficiently large. 910 In the 
case of electromyography, several thousand temporal and 
linear data points are used for the analysis; however, only 
tens of spatial but two-dimensional data points are used 
because the pixel data from the region of interest (ROI) in 
each muscle must be calculated. FFT of EI would generate 
a data structure with a lower sensitivity and larger errors 
than would electromyography. 

The aim of this study was to investigate and compare 
the motor functions and EI of related muscles in young and 
elderly individuals and to analyze the frequency component 
of EI using the maximum entropy method, which has primar- 
ily been used to analyze heart rate variability. 1112 However, 
using this method, frequency components can be calculated 
without the window function, which is appropriate for 
analyses of limited data. A newly designed muscle quality 
evaluation system was proposed in this study. 

Methods 

Subjects 

We enrolled 1 9 elderly individuals aged s65 years (mean age + 
standard deviation, 73+3.2 years; range, 65-77 years; 14 males, 
5 females) with no history of neuromuscular disease and nor- 
mal participation in the activities of daily living as well as 



19 young university students (mean age, 22+1 .5 years; range, 
21-26 years; 14 males, 5 females) (Table 1). Approval for the 
study was obtained from the Saitama Prefectural University 
Ethics Committee. Only individuals who provided informed 
consent for participation were included. 

Motor function measurement 

Motor function measurements were collected in accordance 
with the "Appraisal standards for functional measures in 
independent elderly people" set by the Ministry of Health, 
Labor and Welfare, Japan. The method was demonstrated to 
the subjects before beginning the procedure, and the measure- 
ment process was begun only after the subjects adequately 
understood the procedure and after sufficient recovery peri- 
ods were provided after the preliminary training and between 
the first and second trials. Four areas of motor function were 
measured as follows: 

1 . Normal walking (5 m): the subject was asked to walk for 
a total distance >8 m. The time interval between when a 
foot first crossed the 3 m mark until the foot crossed the 
end of the measurement area was recorded. The standard 
instruction given to all subjects was to "Walk normally 
in a straight line." 

2. Fastest walking (5 m): the subject was again asked to 
walk for a total distance >8 m; however, the standard 
instruction given to all subjects was "Walk as quickly as 
possible in a straight line." 

3. Timed up and go (TUG) test: the subject began the test 
while seated on a chair with his/her back pressed against 
the back support; the time required to stand, walk to 
a marker placed 3 m away, turn the marker over, and 
return to a seated position was measured. The standard 
instruction given to all subjects was "Return to the seat 
as quickly as possible." 

4. Muscle strength: the maximum isometric knee extension 
torque was measured in the dominant leg. The subject was 
seated in a chair with his/her hip and knee in the initial 



Table I Subject characteristics 



Variables 


Elderly (n=l9) 


Young (n=l9) 




male (n=l 4), 


male (n=l4), 




female (n=5) 


female (n=5) 




Mean ± SD 


Mean ± SD 


Age (years) 


73.0±3.2 


22.0±l.5 


Height (cm) 


I60.0±7.7 


I68.9±0.0 


Body mass (kg) 


59.7±9.6 


63.2±9.2 


BMI (kg/m 2 ) 


23.2±2.6 


22.1 ±2.4 



Abbreviations: BMI, body mass index; SD, standard deviation. 
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position of 90° flexion. A dynamometer for isometric 
contraction (|i Tas F-l, operation range 0-980 N; 
Anima, Tokyo, Japan) was used to measure the force. The 
subject's ankle was attached to a belt that was fixed to 
both the dynamometer and the chair. The subject exerted 
maximum effort in an attempt to extend his/her knee 
for >3 seconds, during which the maximum exerted value 
was measured. The length between the lateral knee joint 
fissure gap and ankle attachment was measured. The knee 
extension torque was calculated as the strength multiplied 
by the length. 

High intra-class correlation coefficients (ICCs) more 
than 0.95 have been reported for similar test with the normal 
walking, fastest walking, and TUG tests," 14 and for similar 
tests for muscle strength. 1 "' 

Each motor function area was measured twice per subject, 
and the results were averaged. 

Ultrasound image recording 

The cross-sectional echo images were recorded using a 
B-mode ultrasound device (Mysono U6; Samsung Medison, 
Seoul, Korea). A linear probe (5-12 MHz) was set to 8.8 
MHz and placed on the skin perpendicular to the skin surface 
(Figure 1). All system settings and parameters remained con- 
stant for all subjects. The time-gain compensations according 
to the tissue depths were set to neutral. The fully relaxed 
subject was seated in a chair with his/her hip and knee at 
90° flexion. The muscle strength measurement position was 
the midpoint between the lateral epicondyle of the femur and 
the ipsilateral greater trochanter. The probe was positioned 
perpendicularly to the longitudinal quadriceps femoris (QF) 
muscle axis, and a transverse plane scan was performed at a 
depth of 5 cm. The probe head was coated with a transmis- 
sion gel that provided acoustical contact without compressing 
the skin surface. 



around the central areas of RF and VI without multiple 
reflection areas to obtain ROIs of RF and VI. ROI pixel data 
were converted to a comma-separated value (CSV) format 
numeric table (64x64 row-column) using Origin Pro version 
8.6.0 (Origin, Wheeling, IL, USA). EI of each muscle was 
determined using the mean grayscale values from 0-256, 
expressed in ROI values (0: black; 256: white). 

The mean frequency of each column was analyzed 
according to the maximum entropy method using MemCalc 
version 2.0 (Suwa Trust, Tokyo, Japan). 16 The grayscale 
values of each column comprised the superficial-to-deep 
64 rows. The analyzed data comprised a space series of a 
row x column-pixel image of a single instant, not the tem- 
poral series (fluctuated data with the passage of time) gener- 
ally used for this type of frequency analysis. In this study, 
the space interval between neighboring pixels (or grayscale 
values) was 0. 1 mm. In MemCalc, a space series is assumed 
to comprise an underlying variation and fluctuating parts; 
the underlying variation is expressed as the function X uy {s), 
which can be described as a linear combination of sine and 
cosine functions, as follows: 

Np 

XwW = flo+Sk sm(2xf n s)+b n cos(2^/>)] (1) 

where^ is the frequency of the nth component, a n and b n are 
the amplitudes of the «th periodic component, N is the total 
number of components, and a 0 is a constant that indicates 
the mean time series value. The value off is determined by 
peaks in the power spectral density. The estimated value P(f) 
can be expressed as follows: 



P(f)-- 



AsP 




1+ S r mk e*P(-&*JMs) 



(2) 



Ultrasound image analysis 

An electric caliper was used to measure rectus femoris (RF) 
and vastus intermedius (VI) thicknesses on the recorded 
ultrasound images. QF thickness was defined as the sum 
of RF and VI thicknesses. If a fixed scanning depth of 
5 cm was insufficient to display the entire muscle, the depth 
was widened to 5.2 cm for thickness measurements only 
(Figure 1). 

The ultrasound images were exported in bitmap format to 
a personal computer to measure EI and allow more detailed 
analyses. Adobe Photoshop CS6 version 13.0 (Adobe Systems, 
San Jose, CA, USA) was used to extract 64x64-pixel ROIs 



where P is the prediction error filter output power of 
the order m and y , is the corresponding filter coefficient 
m =0, 1, 2, M (M = optimum filter order). P and y 
are determined by Yule-Walker equations, using Burg's 
algorithm. The details of these equations and the optimum 
setting values have been described in previous references, in 
which this method was used to analyze temporal heart rate 
variability data. 1617 The calculated mean frequencies of the 
64 columns were averaged as the mean frequency of ROI 
(MFROI; Figure 1). 

High ICCs more than 0.90 have been reported for the 
muscle thickness and the EI. 1S 
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Figure I Example MFROI calculations from ultrasound images. 

Notes: The mean frequency of each column in a 64x64-pixel ROI (upper half) is calculated and the mean ROI frequencies are then averaged to obtain each MFROI (lower 
half). Notice that the frequencies are distributed to lower values ("frequency distribution of each column"). 

Abbreviations: El, echo intensity; MF, mean frequency; MFROI, mean frequency of region of interest; RF, rectus femoris muscle; ROI, region of interest; VI, vastus 
intermedius muscle. 



Statistical analyses 

Statistical processing was conducted with SPSS version 17.0 
(SPSS Inc., Chicago, IL, USA). The measured values have 
been reported as means + standard deviations. Statistical sig- 
nificance was defined as a P-value <0.05. The non-parametric 
independent two-sample Mann- Whitney (7-test was used to 
compare the motor functions and ultrasound images of elderly 
and young subjects; the non-parametric dependent two-sample 
Wilcoxon test was used to compare the ultrasound images of 



RF and VI ROIs. In both elderly and young subjects, Pearson 
correlation coefficients and two-tailed significance were used 
to estimate the correlations between the measured data of the 
areas with respect to the motor function, EI, and MFROI. 

Results 

Motor function comparison 

Compared with young subjects, elderly subjects had reduced 
muscle strength (P<0.01) and slower 5 m fastest walking 
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times OP<0.001) and TUG test times (P<0.001). There was 
no significant difference in the 5 m normal walking times 
(Table 2). 

Ultrasound image comparison 

Compared with young subjects, elderly subjects had reduced 
RF (P<0.001), VI (P<0.05), and QF thicknesses (P<0.001 ; 
Table 2). EI values of RF and VI were significantly greater 
in elderly subjects than in young subjects (P<0.01 and 
P<0.001, respectively). MFROIs of RF and VI were also 
significantly greater in elderly subjects than in young subjects 
(P<0.01 andP<0.05, respectively). In a comparison of ROI 
Els, RF EI significantly differed from the VI EI in young 
subjects but not elderly subjects (Figure 2A). In a comparison 
of ROI MFROIs, RF MFROI did not significantly differ from 
VI MFROI in both elderly and young subjects (Figure 2B). 

Correlations between the motor 
functions, muscle thicknesses, El, 
and MFROI 

A significantly positive correlation was identified between 
fastest walk and normal walk in both elderly subjects and 
young subjects, and TUG exhibited significant positive cor- 
relations with normal walk and fastest walk in young subjects 
(Table 3). A significantly positive correlation was identified 
between RF thickness and muscle strength in both elderly 
subjects and young subjects, but no significant correlations 
were found between the other motor function and ultrasound 
image values, and QF thicknesses exhibited significant posi- 
tive correlations with RF and VI thicknesses in both elderly 



Table 2 Subjects' motor functions and ultrasound image analysis 
results 



Variables 


Elderly 
Mean ± SD 


Young 
Mean ± SD 


P-value 


Normal walk (sec) 


3.4S±0.37 


3.40±0.55 


NS* 


Fastest walk (sec) 


2.64±0.26 


1 .95±0.38 


<0.00l 


TUG test (sec) 


5.38±0.67 


4.08±0.55 


<0.00l 


Muscle strength (Nm) 


72.92±3I.3I 


I05.36±29.02 


<0.0I 


RF muscle thickness (mm) 


I4.30±2.74 


2I.65±3.54 


<0.00l 


VI muscle thickness (mm) 


14.31 ±2.45 


1 6. 1 5+2.34 


<0.05 


QF muscle thickness (mm) 


28.61 ±3.56 


37.80±4.30 


<0.00l 


RF El (AU) 


25.29±l 1.27 


I6.49±5.22 


<0.0I 


VI El (AU) 


2I.88±9.63 


9.95±4.78 


<0.00l 


RF MFROI (mm" 1 ) 


0.79±0.I3 


0.60±0.22 


<0.0I 


VI MFROI (mm" 1 ) 


0.85±0.I3 


0.70±0.20 


<0.05 



Note: *Not significant (P>0.05). 

Abbreviations: sec, seconds; AU, arbitrary unit; NS, not significant; El, echo intensity; 
MFROI, mean frequency of region of interest; QF, quadriceps femoris; RF, rectus 
femoris; SD, standard deviation; TUG, timed up and go; VI, vastus intermedius. 



subjects and young subjects. VI Els exhibited significant 
negative correlations with VI and QF thicknesses in both 
elderly subjects and young subjects. RF MFROIs exhibited 
a significant negative correlation with VI and QF thicknesses 
and a positive correlation with VI EI in elderly subjects. 

Discussion 

The most novel proposal of this study was that of a method 
designed specifically to investigate muscle tissue quality 
using EI frequency analysis. Using this method, this study 
yielded three important findings. First, despite a lack of 
significant difference in normal walking, the healthy elderly 
subjects had significantly larger MFROI values than the 
young subjects. Secondly, the MFROIs did not appear to 
differ greatly between the superficial RF and the deeper VI 
of ROI relative to the Els because RF EI was greater than VI 
EI in the elderly subjects but not significantly different in the 
young subjects; however, RF MFROIs did not significantly 
differ between the elderly and the young subjects. Thirdly, the 
relationships between the MFROIs and muscle thicknesses 
differed in the elderly and the young subjects because RF 
MFROIs exhibited a significant negative correlation with 
VI and RF thicknesses and VI EI in the elderly subjects; 
however, no significant correlations were observed with any 
analyzed echo image values in the young subjects. 

El value detection 

This study investigated the skeletal muscle quality in healthy 
elderly and young subjects with grayscale analysis in a more 
detailed method compared to that used in previous studies. 
An identical ROI for each RF and VI was used to analyze 
EI in all subjects. The thinnest measured muscle was suf- 
ficiently thicker than the 64-pixel columns. The pixel data, 
which were collected near the surrounding muscle fascia, 
were not included for any of the subjects. 

Li et al evaluated the Els of 40x40-pixel ROIs around 
the forearm flexor muscles and reported that the superficial 
tissue EI values were significantly greater than those of 
deeper tissues in young subjects but that no such difference 
was observed in elderly subjects. 19 Similarly, in this study, 
the RF EI value, which is located more superficially than VI 
EI, was significantly greater than VI EI in the young sub- 
jects but not in the elderly subjects (Figure 1 ; Figure 2). The 
decrease in EI according to the tissue depth in the ultrasound 
images did not correspond with EI attenuation because the 
ultrasound device had been contrived to improve the image 
quality regardless of the tissue depth. EI is influenced by the 
ultrasound device settings, meaning that the grayscale values 
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A A comparison of Els B A comparison of MFROIs 




Elderly Young Elderly Young 



Figure 2 Comparison of the calculated El and MFROI values in the different muscles. 

Notes: RF El is significantly larger than VI El in young subjects (A); however, RF MFROI does not significantly differ from VI MFROI (B). **P<0.0l, NS, not significant 
(P>0.05). 

Abbreviations: Els, echo intensities; MFROI, mean frequency of region of interest; RF, rectus femoris muscle; VI, vastus intermedius muscle. 



cannot be generalized by simple value averaging. Therefore, 
Pillen et al attempted to normalize EI values to improve the 
reproducibility. 20 

Validation of the calculated MFROI 

The probe scanning frequency is closely related to the 
echo image resolution and echo attenuation. The average 
sound velocity through human soft tissues is 1,540 m/s, 
which is the calibrated velocity for most ultrasound 
devices. 21 The axial (ultrasound image column direction) 
EI frequency for the scanning of ultrasound images in 
this study was theoretically estimated as a maximum of 
2.86 mnr 1 with a probe frequency of 8.8 MHz (wavelength: 
l/( 1 ,540 m/s /8,800,000 Hz) =0.000 1 75 m =0. 1 75 mm; reso- 
lution: 1/0.175 mm =5.71 mm" 1 ; maximum EI frequency: 
5.71 mm -1 x0.5 (Nyquist theorem) =2.86 mnr 1 ). The real 
calculated EI frequencies correlated well with the theoreti- 
cal values because most of the frequencies had distributions 
<2.86 mnr 1 in all subjects (Figure 1; Table 2). The space 
interval of 0. 1 mm indicates a space sampling frequency of 
10 mm -1 (1/0.1 mm =10 mnr 1 ). Therefore, the frequency 
was calculated from 0-5 mnr 1 (10 mnr 1 x0.5 [Nyquist 
theorem] =5 mnr 1 ) in this study ("Frequency distribution 
of each column" in Figure 1). Only the column direction 
frequencies were analyzed and averaged because in ultra- 
sound imaging, the row direction resolution is inferior to 
the column direction resolution. 

As a pilot test, the correlation coefficients were calcu- 
lated using two images taken from the left and right halves 
of the 64x64-pixel ROIs and confirmed a high correlation 
(not shown in this study). More validation of this MFROI 



calculation method would also be needed to confirm the 
reproducibility as well as showing any measurement bias 
for further study. 

Estimated advantage of MFROI 
calculation 

RF and VI MFROIs in elderly subjects were significantly 
greater than those in young subjects. However, the calcu- 
lated MFROIs did not differ between RF and VI Els in both 
elderly and young subjects. MFROI might also be more 
depth-independent than EI. The mean EI frequency calcula- 
tion could be a general method that is not largely influenced 
by the depth of the tissue in an ultrasound image. 

In ultrasound images, the anterior thigh muscle usually 
appears as predominantly hypoechoic, although interspersed 
hyperechoic septa can be identified. 22 The calculated MFROI 
would therefore reflect repetitions of the muscle fibers or 
fascicles and the septa. The large value of MFROI in an 
elderly individual would be related to increased levels of 
diffused adipose tissues in the muscle fibers or fascicles; 
this has been confirmed using computerized tomography and 
magnetic resonance imaging. 12 - 23 - 24 The smaller muscle fibers 
or fascicles consequent to age-related muscular dystrophy 
might increase MFROI in elderly individuals. 

The MFROI and EI could be caused by more age-related 
factors than the muscle strengths or muscle thicknesses 
because the muscle strengths exhibited no significant cor- 
relations with the MFROI or EI, although they exhibited 
significant correlations with the muscle thicknesses in both 
elderly and young individuals whose ages were narrowly 
ranged in each group. 
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Notes: *P<0.05, **P<0.0 1 . 

Abbreviations: El, echo intensity; FW, fastest walk; MF, mean frequency of region of interest; MS, muscle strength; MT, muscle thickness; NW, normal walk; QF, quadriceps 
femoris muscle; RF, rectus femoris muscle; TUG, timed up and go test; VI, vastus intermedius muscle. 



We conclude that the ultrasound imaging EI increased 
in elderly subjects as a result of adipose and connective 
tissue accumulation that induced changes to frequencies 
higher than around 0.65 cycles/mm (the mean value of RF 
and VI MFROI in young subjects was 0.65 mm -1 ), although 
elderly subjects did not differ from healthy young subjects 
when walking at a normal speed. The negative correlation 
between QF thickness and RF MFROI suggests that age- 
related muscle mass decreases along with adipose and/or 
connective tissue accumulation, a correlation that was not 
observed in young subjects. 

Based on these results, the main advantage of ultrasound 
image frequency analysis is apparent: the change in EI fre- 
quencies per mm could be directly related to the anatomical 
tissue structure, with only small influences from the scanning 
conditions of an ultrasound device. 

These newly suggested methods could be useful for an 
easily accessible and safe method to evaluate muscle quality 
with aging that may offer more detailed information than 
only EI, and further investigations of muscle quality and 
applications for the early prevention of age-related motor 
functional declines. 
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